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Abstract. We study the interaction problem of a linear polymer chain, floating 
in fractal containers that belong to the three-dimensional Sierpinski gasket (3D SG) 
family of fractals, with a surface-adsorbed linear polymer chain. Each member of the 
3D SG fractal family has a fractal impenetrable 2D adsorbing surface, which appears 
to be 2D SG fractal. The two-polymer system is modelled by two mutually crossing 
self-avoiding walks. By applying the Monte Carlo Renormalization Group (MCRG) 
method, we calculate the critical exponents tp, associated with the number of contacts 
of the 3D SG floating polymer chain, and the 2D SG adsorbed polymer chain, for 
a sequence of SG fractals with 2 < b < 40. Besides, we propose the codimension 
additivity (CA) argument formula for <p, and compare its predictions with our reliable 
set of the MCRG data. We find that ip monotonically decreases with increasing 6, 
that is, with increase of the container fractal dimension. Finally, we discuss the 
relations between different contact exponents, and analyze their possible behaviour 
in the fractal-to-Euclidean crossover region b — > oo. 
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1. Introduction 

The self-avoiding walk (SAW) is a random walk that must not contain self-intersections. 
This kind of random walk, placed on a lattice, has been widely used as a model of a 
linear polymer in a good solvent [1] . Even though an isolated polymer chain is difficult to 
observe experimentally, numerous studies of the single-chain statistics have been upheld 
as a requisite step towards understanding the statistics of collection-chain systems [2]. 
A reasonable extension of a single polymer concept is a model of two interacting 
linear polymers (two interacting SAWs), which has been a subject of extensive studies 
because of its theoretical and practical interest. To investigate the critical properties 
of the two-chain (or many-chain) systems various theoretical techniques have been 
applied, including the field theoretical approach [3], Monte Carlo simulations [4,5], 
transfer-matrix calculations [6], and renormalization group (RG) methods [7-9]. Among 
other fields, a system of two interacting SAWs, that are mutually avoiding, has been 
successfully applied as a model of diblock copolymers [10, 11], as well as in the studies 
of double-stranded DNA molecules [12-14]. 

In order to study the number of contacts between monomers that belong to different 
polymers, the two-polymer system may be modelled by two mutually crossing self- 
avoiding walks [15,16], that is, by two SAWs whose paths on an underlying lattice can 
cross (intersect) each other. Each crossing between two SAW paths corresponds to a 
contact of two monomers that belong to different polymer chains, and therefore with 
each crossing we may associate the contact energy e c . In analogy with the problem 
of polymer adsorption onto a hyperplane [17], we may assume that with decreasing of 
temperature the number of crossings M increases so that at the critical temperature T c 
it behaves according to the power law 

M ~ N v , (1) 

where N is the total number of monomers in the longer chain, and ip is the crossover 
critical exponent. Below T c the number of contacts becomes proportional to N, whereas 
above T c it is vanishingly small. 

In spite of different approaches to the two-SAW problem, the entire physical picture 
achieved so far, is almost entirely of a phenomenological character and there are only 
few numerical results. For instance, in the case of Euclidean lattices, besides the 
codimension additivity (CA) argument predictions for the contact critical exponents 
(with conjectures <p = 1/2 in d — 2, and ip ~ 1/5 in d — 3), there is only Monte Carlo 
result ip = 0.516 ± 0.005, for a model of two mutually crossing SAWs on the square 
lattice [16]. The above problem has been also studied on fractal lattices. But, in the 
case of fractals it may be noted that the two-dimensional (2D) fractal lattices (embedded 
in the 2D Euclidean space) have been more frequently investigated [15,18,19] than the 
corresponding model in the case of 3D fractal structures [20]. Since the 3D fractals 
definitely may serve as a better description of real systems (porous media, for instance), 
it is recommendable to study the described model in a case of a set of 3D fractals. For 
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this reason, it is desirable to extend the study of two-polymer system on a family of 
fractal lattices whose characteristics approach properties of a 3D Euclidean lattice. 

In this paper we report results of our study for the contact critical exponent <p of two 
polymer chains, that display inter-chain interactions, in the three-dimensional fractal 
lattices that belong to the Sierpinski gasket (SG) family of fractals. The two-polymer 
system is modelled by two SAWs which are allowed to cross each other. We assume 
that the first polymer is a floating chain in the bulk of 3D SG fractal, while the second 
polymer chain is adsorbed onto one of the four surfaces (which is, in fact, 2D SG fractal) 
by which is 3D SG fractal bounded. In section[2]of the paper, we first describe the 3D SG 
fractals for general b. Then, we present the framework of the RG method for studying 
statistics of two mutually crossing SAW chains, taking into account the presence of 
the inter-chain interactions, in a way that should make the method transparent for the 
Monte Carlo calculations of the contact critical exponent ip. In section |3] we present the 
obtained specific values of ip for a sequence of 3D SG fractals, that is, for 2 < b < 40. 
In the same section we propose the phenomenological formulae for ip, based on the CA 
arguments, test their predictions on the obtained MCRG data, and discuss relations 
between different contact exponents. Summary of the obtained results and the relevant 
conclusions are given in section HI 

2. Monte Carlo renormalization group approach 

In this section we are going to apply the Monte Carlo renormalization group (MCRG) 
method to the studied model of interacting polymers on the 3D SG family of fractals. 
These fractals have been studied in numerous papers so far, and consequently we shall 
give here only a requisite brief account of their basic properties. It starts with recalling 
the fact that each member of the 3D SG fractal family is labelled by an integer b > 2 
and can be constructed in stages. At the first stage (r = 1) of the construction there is a 
tetrahedron of base b that contains b(b + 1) (6 + 2)/6 upward oriented unit tetrahedrons. 
The subsequent fractal stages are constructed recursively, so that the complete self- 
similar fractal lattice can be obtained as the result of an infinite iterative process of 
successive (r — > r + 1) enlarging the fractal structure b times, and replacing the smallest 
parts of enlarged structure with the initial structure r = 1 (see for instance, figure 1 
of [21]). Fractal dimension df of the 3D SG fractal is equal to 

df = \n[b(b+ l)(6 + 2)/6]/ln6. (2) 

We assume here that one of the four boundaries of the 3D SG fractal is an impenetrable 
adsorbing surface, which is itself a 2D SG fractal with the fractal dimension 

4 D = ln[6(6 + l)/2]/ln6. (3) 

In the terminology that applies to the SAW, we assign the weight x 3 to a step of 
the SAW in the bulk (3D SG fractal), which represents a floating polymer (we mark 
it by P 3 ), and the weight x 2 to a step of the SAW adsorbed on the surface (2D SG 
fractal), which represents an adsorbed polymer (marked by P2), whose monomers act 
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Figure 1. The structure of the three-dimensional SG fractal, for b = 2, at the first 
stage of construction, with an example of the bulk polymer chain (P3) depicted by 
yellow line and the surface-adsorbed polymer chain (P2) depicted by red line. The 
shaded area represents the adsorbing surface (the two-dimensional SG fractal). The 
two polymers are crosslinked at the two cites, so that each contact contributes the 
weight factors w. The blue bonds, marked by i, describe the interactions between those 
monomers which are nearest neighbors to the crosslinked points. Thus, the depicted 
two-polymer configuration contributes the weight x^x^'uP't 3 in the corresponding RG 
equation (especially, in the equation ([7]) for i — 2 and r = 0). 



as pining cites for the floating polymer. In order to explore interacting effects of two 
SAWs on the underlying fractals, we introduce the two Boltzmann factors w = e _£c//fcsT 
and t = e~ et / fesT , where e c is the energy of two monomers in contact (which occurs at a 
crossing site of SAWs), while e t is the energy associated with two sites which are nearest 
neighbours to a crrosslinked site and which are visited by different SAWs (see figure [Q. 

To describe all possible configurations of the two-chain polymer system, within the 
accepted model, we need to introduce the nine restricted partition functions A^ r \ B^ r \ 
CW ; a { [ ] (i = 1,2,3,4), and b\ t) (i = 1,2), that are depicted in figured The recursive 
nature of the fractal construction implies the following recursion relations for restricted 
partition functions 

= £ a(N A , Nb) A Na B Nb , (4) 

N A ,N B 

B (r+D = b(N A ,N B )A NA B NB , (5) 

N A ,N B 

C (r+1) = £ c ( Nc ) C N C ^ (g) 
N C 

Af +1) = a*(A0 A Na B Nb C Nc n Aj A * fl B N k Bk , i = 1, 2, 3, 4 , (7) 

Af j=l k=l 

Bt +l) =Y.H^)A Na B Nb C Nc X{A^\{B^ b \ 1 = 1,2, (8) 

Af j=l k=l 

where Af denotes the set of numbers A/" = {N A , Nb, Nc, N Al , Na 2 , Aa 3 , N Ai , Nb 1} Nb 2 }, 
and, where we have omitted the superscript (r) on the right-hand side of the above 
relations. The self-similarity of the fractals implies that the sets of coefficients a, b 
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Figure 2. Schematic representation of the nine restricted generating functions used 
in the description of all possible two-SAW configurations, within the r-th stage of the 
3D SG fractal structure. Thus, for example, the A\ represents the configuration when 
both polymers (3D SG floating and 2D SG adsorbed) pass through the same vertices 
of a particular r-th stage of the fractal lattice. The interior details of the r-th stage 
fractal structure, as well as details of the chains, are not shown (for the chains, they 
are manifested by the wiggles of the SAW paths). 

and c, that describe a single polymer configurations (a and b coefficients describe the 
bulk polymer configurations, whereas the c coefficients describe the adsorbed polymer 
configurations), and the sets of coefficients a.; and of the corresponding two-chain 
configurations do not depend on r. Each of the two-chain coefficients (a« and bj) 
represents the number of ways in which the corresponding parts of the two-SAW 
configuration, within the (r + l)-th stage fractal structure, can be comprised of the 
two-SAW configurations within the fractal structures of the next lower order. Because 
of the independence of r, these coefficients can be calculated by studying all two-SAW 
paths within the fractal generator only, that is on the first step of construction r = 1. 

The above set of relations (JI])-®, can be considered as the RG equations for the 
problem under study, with the initial conditions: = 23, = x|, = x 2 , 

— X3X2W 2 , A% = A3 = x 3 X2wt, Af' = x 3 x 2 and = B^ = x|x 2 w 2 , which 
correspond to the unit tetrahedron. On the physical grounds that are reasonable for 
the studied model [22], one can expect that the RG transformations should have three 
relevant fixed points of the general type 

(A*, B*, C*, A*,, A*, A*, At, B{, B*). (9) 

The first fixed point 

(A*, B*, C*, 0, 0, 0, At, 0, 0), (10) 

with A* = (z = 1, 2, 3), A\ = A\(A*, B*, C*) and B* = (z = 1, 2), due to the meaning 
of these quantities (see figure 2), describes segregated phase of two chain polymers that 
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should be expected in the high temperature region T > T c . The second fixed point 

(A*, B*, C*, A*C*, A*C*, A*C*, A*C*, B*C*, B*C*), (11) 

with A* = A*C*(i = 1,2,3,4) and B* = B*C* (i = 1,2), which appears to be a 
tricritical point, describes the state of the two-polymer system that occurs at the critical 
temperature T = T c when segregated and entangled polymer phases become identical. 
Finally, the third fixed point 

(0,0, 0,C*, 0,0, 0,0,0), (12) 

with the nonzero value only for A\ = C*, describes the polymer entangled state, which 
should appear at low temperatures T < T c . 

In what follows we focus our attention on the tricritical fixed point fill I) to calculate 
the contact critical exponent ^32 between P3 and P2 polymer chains. We can observe 
that the above system of RG equations (]4|)-(J8]) can be split into three uncoupled sets of 
RG equations: (JI])-©, ([HD, an d CO) - ®- Moreover, it should be noticed, that for each 
b, the first two sets of RG equations, (JI])-© and ([6]), have only one nontrivial fixed 
point value (A*,B*) [21] and C* [23] respectively, which thereby completely determine 
the coordinates of the tricritical fixed point. Calculation of ^32 starts with solving the 
eigenvalue problem of the RG equations (J4j) — dSJ) , linearized at the tricritical fixed point. 
The related eigenvalue problem can be separated into three parts. The first part of 
the eigenvalue problem, related to the equations gives the eigenvalue \ U3 of the 

end-to-end distance critical exponent u 3 = lnfr/lnA^ of the 3D SG floating polymer, 
while the second part, related to the equation §6§, gives the eigenvalue X U2 of the end- 
to-end distance critical exponent z/ 2 = lnfe/lnA^ of the 2D SG SAW, that represents 
the adsorbed polymer. Finally, the third part of the eigenvalue problem (related to the 
equations (CD)-®) reduces to solving the equation 

' dX - ' 0, (13) 



det 



- A 5ij 



where Xj are elements of the set {A\, A 2 , A 3 , A±, B±, B 2 }, and the asterisk means that 
the derivatives should be taken at the tricritical fixed point. Also, we have used the 
prime symbol as a superscript for the (r + l)-th restricted partition functions and no 
indices for the r-th order partition functions. The largest eigenvalue X V32 of above 
equation determines the contact critical exponent via the formula 

^32 = (14) 

In A„ 3 

Hence, in an exact RG evaluation of (^32 one needs to calculate partial derivatives 
of sums (T4])-([ED, and thereby one should find the coefficients a, b, c, aj, and 6j. The 
latter can be calculated by an exact enumeration of all possible two-SAW configurations 
for each particular b. We have found that this enumeration is feasible only for small 
b. However, for large b the exact enumeration turns out to be a forbidding task. We 
bypass this problem by applying the MCRG method. Within this method, the first step 
would be to locate the tricritical fixed point. To this end, we may observe that the 
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results obtained in [21] provide information for both (A*,B*) and X U3 for a sequence of 
3D SG fractals with 2 < b < 40, whereas the results obtained in [24] give the fixed point 
values of C* for 2D SG fractals, for b in the range 2 < b < 80 . Accordingly, the next 
step in the MCRG method consists of finding without explicit calculation of the 
RG equation coefficients. 

To solve the partial eigenvalue problem ([TBI , so as to learn A^ 32 , we need to find 
the requisite 36 partial derivatives. These derivatives can be related to various averages 
of the numbers N Ai and A?^ of different crossings of the SAWs for various two-SAW 
configurations that correspond to the restricted partition functions and Bf\ For 
instance, starting with (J7J) (in the notation that does not use the superscripts (r + 1) 
and r) and by differentiating it with respect to A\ we get 

^ = Y.N Al a % WA NA B N -C N -A^\{A N ^\{B N k B \ i = l,2,3,4. (15) 
dA i Tr ;=i k=i 

Now, assuming that A[ represents the grand canonical partition functions for the 
ensemble of all possible two-SAW configurations, where each of two SAWs starts and 
leaves fractal generator at two fixed corners, so that the first SAW is a 3D SG floating 
chain, while the other SAW is a 2D SG adsorbed chain. With this concept in mind, we 
can write the corresponding ensemble average 

(N Al ) A , = ± E NaM^) A n -B n *C n c n 4 Aj fl < Bk , .= 1,2,3,4, (16) 

A i M j=l k=l 

which can be directly measured in a Monte Carlo simulation. Combing ffT5l) and ffT6l) 
we can express the requisite partial derivative in terms of the measurable quantity 
8 A' A' 

gf^f^h,, i=l,2,3,4. (17) 

In a similar way we can find the rest derivatives, so that all needed derivatives, calculated 
in the tricritical fixed point, may be related with the corresponding averages 

fBA'\* ( 8A'\* A* 

=«V {ai) ^ = 1,2,3,4; t=l,2, (18) 

(f)R'\* R* / f)R'\* 

BlQ) =*<Wf y=M^ i = UM;M = y. (19) 

Consequently to solve the eigenvalue problem (ITBl) . so as to learn X V32 , we need to find 
the above partial derivatives at the tricritical fixed point. These derivatives are related to 
various averages of the numbers N^., and Ng k , of different two-SAW parts (of the types 
Aj, and B k ) within the corresponding two-SAW configurations (described with A\ or B[). 
Therefore, to calculate the derivatives (TTH]) — (jTU]) at the tricritical fixed point, one needs 
36 averages ((Na 3 )* a >, (N Bk ) A ,, (N Aj )g,, (N Bk )* B ,^, which are all measurable through 
Monte Carlo simulations. The pertinent Monte Carlo technique has been described 
in [24], and we would not like to elaborate on them here. Solving numerically the 
eigenvalue equation (|T3|) we obtain \ V32 , and, finally, using relation (|T4"I) . we calculate 
the contact critical exponent (^32, whose specific values we present in the next section. 
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Table 1. The MCRG (2 < b < 40) results obtained in this work for the contact critical 
exponents ^32 for the 3D SG family of fractals. For the completeness we quote here 
the MCRG values for the RG parameters (A*,B*) [21], and C* [24]. 



2 
3 
4 
5 
G 
7 



.4* 



B* 



C* 











8 

9 

10 
12 
15 
17 
20 
25 
30 
35 
40 



4311±0.0009 
3421±0.0004 
2898±0.0004 
2560±0.0004 
2319±0.0003 
2148±0.0003 
2016±0.0003 
1912±0.0004 
1829±0.0003 
1703±0.0004 
1581±0.0001 
1526±0.0001 
1462±0.0001 
1399±0.0001 
1353±0.0001 
1327±0.0001 
1305±0.0001 



0.0505±0 
0.0245±0 
0.0122±0 
0.0067±0 
0.0038±0 
0.0020±0 
0.0012±0 
0.0007±0 
0.0005±0 
0.0001±0 



0023 
0015 
0020 
0019 
0012 
0018 
0026 
0008 
0023 
0035 



0.61825±0. 
0.55137±0. 
0.50658±0 
0.47455±0. 
0.45091±0. 
0.43240±0. 
0.41780±0, 
0.40574±0. 
0.39586±0, 
0.38037±0. 
0.36396±0, 
0.35593±0, 
0.34681±0. 
0.33602±0 
0.32876±0. 
0.32350±0. 
0.31936±0 



00061 
00044 
00034 
00028 
00024 
00021 
00019 
00017 
00007 
00013 
00011 
00008 
00006 
00008 
00007 
00008 
00006 



0.5440±0.0056 
0.4969±0.0024 
0.4658±0.0006 
0.4451±0.0012 
0.4250±0.0004 
0.4092±0.0008 
0.3963±0.0006 
0.3841±0.0007 
0.3714±0.0005 
0.3514±0.0004 
0.3226±0.0003 
0.3126±0.0003 
0.2956±0.0003 
0.2677±0.0003 
0.2573±0.0002 
0.2183±0.0003 
0.2016±0.0003 



3. Results and discussion 

We have studied the interaction problem of a floating polymer chain confined in the 3D 
SG fractal container, with a surface (2D SG fractal) adsorbed polymer chain to calculate 
the critical exponent ^32 which governs the number of contacts between two chains. 

The main goal of this study is the MCRG evaluation of (^32, for various values of 
b. Furthermore, in the case b = 2, we have calculated the exact value. To find the 
exact forms of RG equations (J4j) — (ISI) , for b = 2 fractal, using the computer facilities 
we have been able to enumerate all coefficients (which are available upon request 
addressed to the author) that appear in these equations. Linearization of obtained 
RG equations around the tricritical fixed point (II ip . in this case determined by the 
values A* = 0.4294, B* = 0.0499 and C* = 0.6180 [25], gives the following eigenvalues: 
\ U3 (b = 2) = 2.7965 and X m2 (b = 2) = 1.7475, whereupon we find the exact value 
V 9 32(^ — 2) = In 1.7475/ In 2.7965 = 0.5428. Here, we notice that the same model is 
studied on the four-simplex lattice, which belongs to the same universality class as the 
6 = 2 3D SG fractal. In this study [22] the value (^32 = 0.5669 is reported, which appears 
as an approximate result (and deviates 4% from our exact finding), since in approach 
applied in [22] some two-SAW configurations are not taken into account (more precisely, 
the configurations described by function A% , in our notation). 

For larger b, that is, for fractals in the range 2 < b < 40, we have applied the 
MCRG method expounded in the previous section. The obtained MCRG values for 
y? 32 , together with the pertaining error bars (determined from statistics of measured 
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Figure 3. Results obtained in this work for the contact critical exponent ^32 for 
the two-polymer system, when the first polymer is a floating 3D SG chain, while the 
other is a 2D SG adsorbed chain. The red triangles represent the MCRG values, while 
green triangles correspond to the CA values calculated from ([22]) . The green circle 
depicts the tree-dimensional Euclidean value ip§ 2 ~ 0. The thin solid lines (red and 
green) represent a simple interpolation of data, while the red dotted lines show the 
two possible scenarios for the fractal-to-Euclidean behaviour of 1^32, and serve as the 
guides to the eye. The error bars, for the MCRG data, are not depicted in the figure, 
since in all cases they lie within of the corresponding symbols (red triangles). 



quantities through the Monte Carlo simulations), are given in table [U The values of 
<£>32 for fractals 2 < b < 30, have been calculated from averages obtained from 10 6 
Monte Carlo simulations performed for each possible two-SAW configuration, while 
y?32 for fractals b = 35 and b = 40 was obtained performing 5 10 5 Monte Carlo 
simulations. The time needed to calculate all requisite averages ((Na)a'i (N b )*,, , 
(Naj)* B ', (NB k ) B i) increases exponentially with the scale parameter b, so that for 5 10 5 
simulations performed on fractals b = 10 and b = 40, at the corresponding tricritical 
fixed points, it was required 10 minutes and 90 hours of the available CPU time, 
respectively (on a PC with the Intel Pentium 4 processor). 

Comparing our MCRG and exact result for b = 2 fractal (0.5440 ± 0.0056 versus 
0.5428) we can see that the MCRG result deviates 0.22% from the exact RG finding, 
which is a very well agreement between the two (Monte Carlo and exact) approaches of 
solving the problem, and thereby provides reliance on accuracy of the MCRG results for 
larger b. In figure [3] we depict our MCRG findings for the critical exponent <y9 32 , of the 
3D SG family of fractals, as a function of 1/b. We observe that the critical exponent (P32, 
in the region 2 < b < 40, is a monotonically decreasing function of b, which means that 
the number of polymer contacts (crossings of the SAW paths) decreases with increasing 
of the fractal dimension of the SG fractals. 

In figure [3j we have also presented the numerical values for ip 32 , that follow from 
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Table 2. Numerical values for the contact critical exponents ^32, ^33 and P22 
evaluated from the proposed phenomenological formulae (based on the CA arguments) 
(|2"2")) - (|2"4"1) . The values for ^3 and v-i have been calculated through the MCRG 
simulations performed in [21] and [24], respectively. 
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0.6169 
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0.5953 


0.6689 


8 


0.6130 


0.7742 


0.3805 


0.5887 


0.6658 


9 


0.6087 


0.7722 


0.3738 


0.5855 


0.6622 


10 


0.6048 


0.7698 


0.3690 


0.5833 


0.6603 


12 


0.5987 


0.7659 


0.3608 


0.5792 


0.6571 


15 


0.5933 


0.7620 


0.3497 


0.5711 


0.6529 


17 


0.5899 


0.7590 


0.3455 


0.5683 


0.6523 


20 


0.5869 


0.7560 


0.3384 


0.5621 


0.6506 


25 


0.5817 


0.7516 


0.3317 


0.5577 


0.6495 


30 


0.5795 


0.7481 


0.3249 


0.5502 


0.6491 


35 


0.5759 


0.7457 


0.3212 


0.5490 


0.6480 


40 


0.5755 


0.7434 


0.3158 


0.5416 


0.6479 



the CA argument predictions adopted for the studied model. Namely, the CA argument 
claims that a codimension d — di P of intersection points is a sum of codimensions d— d\ 
and d — d 2 of intersecting objects [16,26]. Here, d is a dimension of embedded space, d ip 
dimension of intersection points, and, d\ and di dimensions of intersecting objects. In 
the studied problem the embedded space is a 3D SG fractal of fractal dimension d 3 j- D , 
while the intersecting objects are linear polymers with the fractal dimensions d\ = I/V3 
and g?2 = 1/^2, whereupon follows 

d ip =(- + -)- df . (20) 

On the other hand, the total number of intersection points behaves according to the 
power law 

M 32 ~ (R N ) dip ~ N usdi " , (21) 

where (Rn) ~ N U3 is the mean end-to-end distance of the longer polymer chain. From 
the latter scaling form, we find the following phenomenological formula for the contact 
critical exponent 

^2 = (l + V f) ~ ^df , (22) 

whose numerical values, evaluated from the attainable values for u 3 and v 2 on SG fractals, 
are listed in table [2j Comparing the values for <£> 32 predicted by the obtained CA 
formula with our convincing set of MCRG data (see figure [3]), we can perceive that in 
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the first part of the examined region (2 < 6 < 40), the CA formula displays satisfactory 
agreement with the corresponding MCRG values, and, for example, for b = 10 fractal 
the two data (CA and MCRG) have almost identical numerical values. For b < 10, 
the CA formula gives some smaller values for ip 32 than the MCRG method (at most 
9%, for 6 = 2), while for b > 10, the CA predicted values get to be larger than the 
corresponding MCRG findings, and, when b increases, a departure of the two sets of 
data becomes unambiguous (for b = 40, the deviation is 57%). Anyway, from both the 
MCRG method and the CA formula ( 1221) . follow that <y9 32 , in the studied region, is a 
monotonically decreasing function of b. 

Besides of ip 32 , we may define the critical exponent ip%3 which describes the number 
of polymer-polymer contacts, when both polymers are floating chains in the 3D SG 
fractal container (that is, when both polymers are of the type P3), and ip 22 when both 
polymers are adsorbed by 2D SG surface (both of the type P 2 ). For these exponents 
the CA arguments give 

<p& = 2 - u 3 df , (23) 

<p% = 2 - u 2 df . (24) 

A simple comparison of the data, from table [2J reveals that the following relationship 
holds 

V?32 < V?33 < V?22 , (25) 

for each particular b, in the range 2 < b < 40. The first inequality means that 
the number of polymer contacts increases when both polymers being situated in the 
bulk of 3D SG fractal, because of growing fractal dimension of a polymer chain 
when it makes a transition from the 2D SG surface to the bulk of 3D SG fractal 
(dj 3 = I/1S3 > dj 2 = 1/^2) • Ori the other hand, in the case when both polymers 
are adsorbed (situated on 2D SG fractal), the fractal dimension of the embedded space 
decreases d"j D < d 3D , and a chance for two polymers to make a contact increases, 
implying the second inequality in fl25l) . 

The relation between ^32 and ^22, can be tested on the reliable sets of MCRG 
values for ip 32 (calculated in this work) and ip 22 (calculated previously [18] in the range 
2 < b < 100). The values for ip 22 , being always larger than the corresponding CA 
values proposed by (1241) . display a monotonic decreasing from ip 22 {b = 2) = 0.7493 up 
to ip 22 (b = 100) = 0.6735 [18]. Their comparison with ip 32 (from tabled]) shows that 
the relation ip 32 < ip 22 is always satisfied. The relations between ^33 and the other 
two contact critical exponents (^32 and ^22) can be numerically proved only for 6 = 2 
fractal, where the RG value is known ^33(6 = 2) = 0.6635 [20], and from which the CA 
value 0.6516 deviates only 1.8%. To complete the inspection of relation fl25l) and test 
CA formula fl23|) for larger 6, one needs to calculate y3 33 values for 6 > 2, which appears 
to be much complicated task than the calculation of ^32, and may be a topic for a future 
study. 
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Finally, we discuss the possible behaviour of the contact critical exponents in the 
fractal-to- Euclidean crossover region, that is, in the limit b — > oo (d 3 j- D — ► 3, d"j D — ► 2). 
For a SAW situated in Euclidean spaces, the end-to-end distance critical exponent in 
d = 2 takes the exact value vf = 3/4 [27], that was numerically confirmed [28] for 
the triangular lattice (to which the 2D SG fractals approach, when b — ► oo), whereas 
in d = 3, at present, the most accurate evaluates are vf = 0.5874 ± 0.0002 [29] and 
if = 0.58765 ± 0.00020 [30]. Consequently, from the formulas (|22j)-(|21j), we compute 
the Euclidean values: y?f 2 ~ 0, y?f 3 ~ 1/5 and tp 2 2 = 1/2. In the case of 2D SG fractals 
(when both polymers are of the type Pg), using the finite-size scaling arguments it was 
argued [20], that in the limit b — > oo, the asymptotic behaviour of y?22 is described by 
the relation <£> 2 2 — 2 — v-tfPP , which coincides with the CA formula ff24|) . Knowing, 
from the finite-size scaling analysis [31], that v 2 for very large b goes, from below, to the 
Euclidean value 3/4, from the latter asymptotic relation follows that ^22 approaches, 
from above, the Euclidean value (p 22 = 1/2, when b — > 00. 

As regards the obtained MCRG results for <^ 32 , it is hard to say what happens 
beyond b = 40, and specially, to establish whether the critical exponent ^32, in the 
limit b — > 00, goes to the zero Euclidean value, or some non-Euclidean value (see figure 
EJ). To answer this question properly, another method may be needed, for instance, the 
finite-size scaling approach that has been applied in the case of 2D SG fractals [20,31]. 
Here, we can only discuss possible predictions, about the limiting value for <£>3 2 , that 
follow from the corresponding CA formula, for which we need to know the asymptotic 
behaviour of 1/3. Unfortunately, up today, the behaviour of z/3, in the fractal-to- Euclidean 
crossover region b — > 00, exists as an unsolved problem. At this moment, we may say 
that if v-& — > if (which means, see data from table El that 1/3 is a non-monotonic 
function of b), then (122]) gives ^32 — > (p§ 2 ~ 0. On the contrary, if 1/3 tends to some non- 
Euclidean value less than , then <£> 32 will go to some value larger than the Euclidean 
value ip§ 2 . The possibility that z/3, for b — ► 00, tends to some value larger than 
should imply a negative value of 9932 (and consequently, in this case, we cannot exclude 
a situation in which the studied transition turns into first order, with increasing b). 
Similar conclusions may be deduced about the limiting value for the contact critical 
exponent ^33. Particulary, from ( |23l) follows: if u 3 — > then <^ 33 — > ip^ 3 w 1/5, or, 
if V3 approaches the non-Euclidean value, then (^33 also approaches a non-Euclidean 
value. From the exposed analysis we may infer that the behaviour of the contact 
critical exponents, in the fractal-to-Euclidean crossover region, is closely related to the 
corresponding behaviour of the end-to-end distance critical exponents. 

4. Summary 

In this paper we have studied the two interacting linear polymer chains, modelled by 
two mutually crossing SAWs, situated on fractal structures represented by the three- 
dimensional (3D) Sierpinski gasket (SG) family of fractals. We take on that the first 
polymer (P3) is a floating chain in the bulk of 3D SG fractal, while the second polymer 
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chain (P 2 ) is adsorbed onto one of the four boundaries of the 3D SG fractal, which 
appears to be 2D SG fractal. Specifically, we have calculated the contact critical 
exponent ^32, associated with the number of monomer- monomer contacts between 
polymers P3 and P 2 . 

By applying the renormalization group (RG) method, we have calculated the exact 
value of the critical exponent y2 32 , for the first member b = 2 of 3D SG fractal family. 
The specific accomplishment in the course of this work is the calculation of a long 
sequence of values of </? 32 , for 2 < b < 40, obtained by applying the Monte Carlo 
renormalization group (MCRG) method. Our results demonstrate that </? 32 , for the 
studied values of b, monotonically decrease with b, and, in analogy with the behaviour 
of ip 2 2 (that governs the number of contacts between two adsorbed polymers) it seems 
that for b > 40 the critical exponent y? 32 continues decreasing, and in the limit b — > 00 
tends to the Euclidean value ipf 2 ~ 0. Finally, using the codimension additivity (CA) 
arguments we have proposed the phenomenological formulae for the considered contact 
critical exponents, and we have tested their predictions on the obtained sets of data. We 
find that, for fractals labelled by smaller values of 6, the CA proposals give satisfactory 
agreement with existing convincing results. 

On the exposed grounds of the presented investigation we may conclude that the 
set of obtained results of the studied problem has been significantly extended. We 
have demonstrated that the statistics of two crosslinked polymer chains on the family 
of 3D SG fractals can be rewardingly studied by the MCRG method. In particular, 
the MCRG study of the contact critical exponents revealed their interesting behaviour 
as the functions of fractal scaling parameter b, making a step forward to prescribe 
the behaviour of SAW critical exponents in the fractal-to-Euclidean crossover region. 
As a further investigation one may attempt to extend our study for calculating the 
contact critical exponent ip 33 , when both polymers are floating chains in the bulk of 3D 
SG fractal. In addition, to make the studied model more realistic (but more difficult 
to study), it can be supplemented by additional interactions. For instance, one can 
introduce the interactions between the bulk floating chain and the sites of adsorbing 
surface (to promote the adsorbed phase of floating chain), as well as the intra-chain 
interactions (to promote a collapsing phase of the bulk floating chain), and then, in the 
space of interaction parameters the phase diagram can be investigated, together with 
new contact critical exponents at appropriate phase transition fixed points. 
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